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Abstract 

In this paper, Optimal Homotopy Perturbation Method (OHPM) is employed to 
determine an analytic approximate solutions for nonlinear MHD Jeffery-Hamel 
flow and heat transfer problem. The Navier-Stokes equations, taking into ac¬ 
count Maxwell’s electromagnetism and heat transfer lead to two nonlinear or¬ 
dinary differential equations. The obtained results by means of OHPM show a 
very good agreement in comparison with the numerical results and with Homo¬ 
topy Perturbation Method (HPM). 

Keywords: optimal homotopy perturbation method, Jeffery-Hamel, nonlinear 
ordinary differential equations. 


1. Introduction 


5 


The incompressible fluid flow with heat transfer is one of the most applica¬ 
ble cases in various flelds of engineering due to it industrial applications. The 
problem of a viscous fluid between two nonparallel walls meeting at a vertex and 
with a source of sink at the vertex was pioneered by Jeffery [l|, and Hamel [^. 
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Later, the Jeffery-Hamel problem have been studied by several researchers and 
discussed in many textbooks and articles. A stationary problem with a finite 
number of ’’outlets” to infinity in the form of infinite sectors is considered by 
Rivkind and Solonnikov . The problem of steady viscous flow in a convergent 
10 channel is analyzed analytically and numerically for small, moderately large and 
asymptotically large Reynolds numbers over entire range of allowed convergence 
angles by Akulenko et al. Q|. The MHD Jeffery-Hamel problem is solved by 
Makinde and Mhone Q using a special type of Hermite-Pade approximation 
semi-numerical approach and by Esmaili et al. Q by applying Adomian decom- 
15 position method. The classical Jeffery-Hamel flow problem is solved by Ganji 
et al. by means of the variational iteration method and homotopy perturba¬ 
tion method, and by Joneidi et al. [sj, by differential transformation method, 
Homotopy Perturbation Method and Homotopy Analysis Method. The classi¬ 
cal Jeffery-Hamel problem was extended in Q to include the effects of external 
20 magnetic field in conducted fluid. Optimal homotopy asymptotic method is 


n 


applied by Marinca and Heri§anu m and by Esmaeilpour and Ganji 


i|n|. 


The 


effect of ma gne tic field and nanoparticle on the Jeffery-Hamel flow are studied 


12| and Numerical treatment using stochastic algorithms is used by 
Raja and Samar Q. 

In general, the problems as Jeffery-Hamel flows and other fluid mechanics 
problems are inherently nonlinear. Excepting a limited number of these prob¬ 
lems most do not have analytical solution. The aim of this paper is to propose 
an accurate approach to the MHD Jeffery-Hamel flow with heat transfer prob¬ 
lem using an analytical technique, namely OHPM [3, Q: Q- Our approach 
does not require a small or large parameter in the governing equations, is based 
on the construction and determination of some auxiliary functions combined 
with a convenient way to optimally control the convergence of the solution. 

2. Problem statement and governing equations 


We consider a system of cylindrical coordinates with a steady flow of an 
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incompressible conducting viscous fluid from a source or sink at channel walls 
lying in planes, with angle 2a, taking into account the effect of electromagnetic 
induction, as shown in Fig. 1, and the heat transfer. 


The continuity equation, the Navier-Stokes e quat ions and energy equation 
in cylindrical coordinates, can be written as [l^ . |l9| . H, 0: 
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where p is the fluid density, P is the pressure, v is the kinematic viscosity, T 
is the temperature, k is the thermal conductivity, Cp is the specific heat at con¬ 
stant pressure, a is the electrical conductivity, Bq is the induced magnetic field 
and the stress components are dehned as 


dUr 2 

= 2— - -dtvu 

( 5 ) 
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By considering the velocity field is only along radial direction i.e. = Q and 
substituting Eqs. ©-([T]) into Eqs. m and ©, the continuity, Navier-Stokes 
and energy equations become: 
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The relevant boundary conditions, due to the symmetry assumption at the 
channel centerline are as follows: 


dUr _ dT _ A 

dip dip ’ 


Ur 


Uc 

r 


at ip = 0 


60 and at the plates making body of channel: 


Ur = 0, T = ^ at if = a 


( 11 ) 


( 12 ) 


where Uc and Tc are the centerline rate of movement and the constant wall 
temperature, respectively. 

From the continuity equation ([5]), one can get 


rur = f{ip) 

65 where f[ip) is an arbitrary function of ip only. 
By integrating Eq. (US it holds that 


(13) 


P{r, ip) = ^f{ip) + P9{r) (14) 

in which g(r) is an arbitrary function of r only. 

Now, defining the dimensionless parameters: 

»? = f. Piv) = ^, (15) 

where Tc is the ambient temperature, and substituting these into Eqs. o and 
([9|) and then eliminating the pressure term, one can put: 


F'" + 2aReFF' + (4 - H)a'^F' = 0 


(16) 


4 


e" + 2a (2a + RePrF) 9 + pPr [{H + 4a2) + F'^ 

75 subject to the boundary conditions 


= 0 


(17) 


F(0) = 1, F'(0) = 0, F(1)=0 


(18) 


0(1) =0, 0'(O) = O 


(19) 


where Re = is the Reynolds number, F[ = yj is the Hartmann number, 
Pr = 1^, f3 = ^ and prime denotes derivative with respect to rj. 


3. Basic ideas of optimal homotopy perturbation method 

To explain the ideas of the optimal homotopy perturbation method, consider 
the non-linear differential equation 


L [u, u\ u", u", Tj] + giif) + N [m, u , u", u", rj] = 0 
that is subject to the initial / boundary condition: 


( 20 ) 


( 21 ) 


where L is a linear operator, g is a known function, N a nonlinear operator, B 
85 is a boundary operator and T is the boundary of the domain of interest 

Q, 


^ 0 . 


We construct the homotopy 


22 | 


R {u,p) = L {u, u', u", u'", rj) -b g{r]) + pN {u, u', u", u'", rj) = 0 


( 22 ) 


for Eq. (EOl), where p is the homotopy parameter, p G [0, 1]. From Eq. (I22[) 
one gets: 


R (u, 0) = L {u, u', u", rt"', g) -b g{g) = 0 
R (ii, 1) = L (u, u', u", u"', g) -b g{g) -b N {u, u', u”, u'", g) =0 


(23) 
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Assuming that the approximate analytical solution of the second-order can 
be expressed in the form 


u{ri) = uq + pui + p‘^U 2 (24) 

and expanding the nonlinear operator N in series, with respect to the parameter 
95 p, we have: 


N {u,u',u",u'",ri) = N (uo,Uo,Mo,'Uq",? 7 ) +p[uiNu (uq, Mq, Mq, Wo",+ 
+u[Nu' (uo, Uq, Uq, Uq, p) + u'{Nu" (uo, Uq, Uq, Uq , p) + (25) 

+u'-['Nu"' (uo, Uq, Uq, Uq", p)] + p^ (u2Nu + u’^Nu' + ...) 

where Fu = By introducing a number of unknown auxiliary functions 

Hi{p,Ck), i = 0, 1 , 2, ... that depend on the variable p and some parameters 
Cfe, k = 1, 2, ..., s, we can construct a new homotopy: 

100 

n iu,p) = L {u, u', u", u'", p) + g{p) + pHoip, Ck)N {uq, u'q, u((, p) + 
p^ [Hi{p,Ck)uiNu (uo) + H2{p,Ck)u[Nfi' (uq) + H3{p,Ck)u'{Na» (uq) + 
+H 4 {p, Ck)u'('Nu"' (uo)] -I- [H^ip, Ck)u2Na (uq) + He{p, Ck)u2Nu' (uq) + ...] 
Equating the coefficients of like powers of p, yields the linear equations: 

L [uo, Uq, Uq, Uq, p] + g{p) = 0, B = 0 (27) 

L{ui) + HQ{p,Ck)N {uq,u'q,u'i^,u'^',p) = 0, b(ui,^^=0 (28) 

L{u 2) + Hi{p,Ck)uiNu {uq,Uq,Uq,Uq',p) + H2ip,Ck)u[Nu> (uq) + 

+H3{p, Ck)u'{Nu» (uq) + Hiip, Ck)u'{'Nu"' (uq) =0, B (^U2, = 0 

105 The functions Hi{p,Ck), i = 0, 1, 2, ... are not unique and can be chosen 
such that the products Hi ■ UjN^ and ujN^ are of the same form. In this way, 
a maximum of only two iterations are required to achieve accurate solutions. 


(26) 
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The unknown parameters Cfc,A: = l, 2, s which appear in the functions 
Hi{r],Ck) can be determined optimally by means of the least-square method, 
110 collocation method, the weighted residuals, the Galerkin method, and so on. 

In this way the solution of Eq. (1^ subject to the initial / boundary con¬ 
dition (ED) can be readily determined. It follows that the basic ideas of our 
procedure are the construction of a new homotopy (1261) , the auxiliary functions 
Hi with parameters Ck that can be determined optimally leading to the conclu- 
115 sion that the convergence of the approximate solutions can be easily controlled. 

4. Application of OHPM to the MHD Jeffery-Hamel flow and heat 
transfer problem 

Let us present the approximate analytic expressions of /(??) and 0{r]) from 
Eqs. (ITBl) - (froi) by means of OHPM. 

120 For Eqs. (HU) and dni), the linear operator is chosen as L {F) = F"', while 
the nonlinear operator is defined as N (F) = 2aFF' -|- (4 — H)a^F', g{r]) = 0. 
The initial approximation Fq is obtained from equation ED 


E"' = 0, Eo(0) = l, F^(0)=0, Eo(l) = 0 (30) 

The solution of Eq. EQl) is hence 

FQ{'q) = l-rf. (31) 

125 On the other hand, from Eq. ED, one obtains 

Nf{F) = 2aReF', Np'{F) = 2aReF F {4, - H)a^ . (32) 

By substituting Eq. (IHTl) into the nonlinear operator N and into Eq. (I32p 
one retrieves: 


A(Eo) = 2V_2(A + H)77, Nf (Fo) =-2Ar), 
Nf' {FQ) = -Ar]^ + A + B 
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where A = 2aRe, B = (4 — H)a^. 

130 Eq. (E51) becomes 

Fl" + Hoiv, Ck) [2V - 2(A + B)r]] = 0, Fi(0) = F{(0) = Fi(l) = 0 (34) 

We choose HQ{r], Ck) = —60Ci where Ci is an unknown parameter and from 
Eq. (l34l) we obtain 


Fi{r]) = 2ACir]^ - 5{A + B)Cir]^ + {3A + 5B)Cif]^ 

Eq. (EH)) can be written in the form 

Ff + 7Ji(ry, Ck)i-2Ar])Fi + H2{v, Ck)i-Ar,^ + kl + S)F( = 0 
F2(0) = F2'(0)=F2(1)=0 

In this case we choose 
Hi{ri, Ck) = + C^rj + C 4 + —^ , H 2 {tj, Ck) = — + — 

135 such that the solution of Eq. (|36|) is given by 


(35) 


(36) 


F2{r]) 


AC 1 C 2 .^11 I 2 ACiC 2 — 5{A-\-B)OiC2 .^10 
495 ' 720 ' 


I 2ylCiC4-5(A+B)CiC3+5A^CiC6 

504 


(3A+5B)CiC2-5(A+B)CiC4+2ACiC5-10(A^+AB)CiC6 8 , 

336 9 “t 

(3A+5B)CiC3-5(A+B)CiC5-5(A=+AB)CiC6+2ACiC7 7 

210 ^ 


(3A+5B)CiC4 + (13A2+25AB+10B^)CiC6-5(A+B)CiC7 6 , 

120 ^ 

(3A+5B)CiC5 ^5 (3A^+8AB+5B^)CiC6 + (3A+5B)CiC7 , ^^^2 

60 I 24 'I '/ 


(37) 


where 

^ = -^1^2 [(i, - TI4 + TT 2 ) ^ + (3M - Th) 5] - 
-^1^3 [(^-^ + ^)^+(^-^)S]- 
-CiC4[(^-3i3 + ^)A+(i-^)F]-CiC5(^ + |f)- 
-CiCe[{^,-^s-m-ro)^"-im + m + h)^B-lB^]-c,Cr{^^ + §) 
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For p = 1 into Eq. (1^ . we obtain the second-order approximate solution, 
using and Eqs. (I5T]) . (1551) and (1571) : 


140 


2ACiC3-5(^+-B)CiC2 ^10 ^ 


720 


Fiv) = Foil) + Flip) + F 2 iri) = + 

, 2AC^Ci-0(A+B)C^C3+bA'^C^C6 , 

504 'I 

I (0A+5B)CiC2-5(A+B)CiC4.+2ACiC5-10(A'^+AB)CiCa 8 , 

I (0A+5B)CiC3-5(A+B)CiC^-b(A'^+AB)CiC6+2ACiCT 7 , 

^ (3A-|-5a)CiC4-l-(13A^-|-257tB+10B^)CiC6-5(A-|-B)CiC7 ^6 _|_ ^2ACi 

- (37t=^-^87^B-^5a")Cl^C6-K3A-^5B)ClC7 5B)Ci p* + [{3A + 5B) + M - 1] p"^ 

Now, we present the approximate analytic solution for Eqs. (HZl) and (USD- 
The linear and nonlinear operators and the function g are, respectively 


I (3A-t-5i3) ^ f-t ’ 
-60- 


FiB) = Q\ 5(^) = -l, 

N{Q) = 1 -h 2aRePrFe + pPr (1 -h 4a'^)F^ + 

The Eq. (l27l) becomes 

0''-l = O, 0o(l)=O, d^(0)=0 

Eq. (1551) has the solution 


eo{p) = ^ii-p^). 


(39) 


(40) 


(41) 


From Eq. (1591) it follows that 

Neie) = 4a2 -h 2aRePrF, Ns>i9) = 0. (42) 

By substituting Eq. (15T1) into Eqs. (1551) and (155t . one gets respectively 


N{eo) = C-2Dp^ + Ep^, Ne{9o)=L + Kp^ ( 43 ) 

where 

C = 1 + 2 q;^ -I- aRePr + A^a^Pr F PrFip 
D = + 2aRePr + 2BPr(2a^ — 1) -|- PrHB 

( 44 ) 

E = aRePr + AfBa^Pr + PrH/B 
L = Aa^ + 2aRePr, K = —2aRePr 


(38) 
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150 


Eq. (E51) can put 


e'i + ha{ii,Cs){C -207^^^+ Et^^) = Q, 0i(1) = 0;(O) = O . (45) 

Choosing ho{r],Cs) = —SOCs into Eq. (HSl) . one obtains 


0i(7?) = C8 [15^(772-l)-5i4(r,4_l)+£;(ry6_i)] . (46) 

Eq. (EH)) can be written in the form 


6 '^ + hi(7?, Ck)iL + Kif)ei = 0, 02(1) = 0^(0) = 0 


and therefore it is natural to choose the auxiliary function hi as 


(47) 


155 


hiiv, Ck) = ^ (Cg + Ciory + + Ci 2 V^ + Ci^rj^) 

From Eq. izi), it can shown that 

02(7?) = -L(15C -5D + E) [iC9(r?2 - 1) + 1^10(773 - 1)] + 

+ [ISLCCg - (15C - 5£i + E){KC9 + LCn)] 

+ [IbLCCio - (ISC' - 5£> + E)(KCio + LC 12 )] 

+ [il5CK - 5DK)C9 + 15LCCii- 
-(15C - 5i4 + E){KCii + LC 13 )] + 

+ [{15CK - 5DL)Cio + 15 LCCi 2 - K{15C - 5D + E)Ci 2 ] (48) 

+ [{LE - 5DK)Cg + il5CK - 5DL)Cii + IbLCCis- 
-K{15C - 5D + E)Ci 3 ]^ + 

+ [{LE - 5DK)Cio + {15CK - 5DL)Ci2] ^ + 

+ [EKCg + {LE - bDK)Cii + {IbCK - 5£>L)C'i3] + 

+ [EKCii + {LE - bDK)Ci 3 ] + EKCi 2 ^^ + 

The second-order approximate solution of Eqs. C2D and (HU is 

0(7?) = 0o(?7) + 0i(7?)02(7?) (49) 

where 0o, 0i and 02 are given by Eqs. (|iT|) . m and (H51) respectively. 
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5. Numerical results 


160 In order to show the efficiency and accuracy of the OHPM , we consider some 
cases for different values of the parameters a and H. In all cases we consider 
Re = 50, Pr = l, 13 = 3.492161428 • IQ-^^ 

Case 5.1 Consider ^ and H = 0. By means of the least-square method, 
the values of the parameters C^, f = 1, 2, ..., 13 are 

Cl = -0.000025737857, C 2 = 128165.4247848388, C 3 = -360860.8730122449, 
C 4 = 315974.73981422884, C 5 = -20823.768289134576, 

Ce = -319.2089575339067, C 7 = -62701.61063351235, 

Cs = -8.904447449602-10■^^ Cg = -2.9799239131772537-10■^^ 

Cio = 0.441808726864- lO'^^ Cn = -1.808778662628 - lO’^^ 

C 12 = 0.111699001609-10-1^ C 13 = -5.293450955214 - lO’^^ 

One get approximate solutions from Eqs. (|38ll and (|4^ respectively: 


F{t]) = 1 - 2.3104494668772 + 2.4868857696774-^ 0.3531718162775- 
- 3.4153413394776 -t 1.7515474936?7^ -f I. 203 I 8 O 5 O 64778 - (50) 

-I. 6209 O 6688 O 776 -t 0.6391440832?7i° - 0.08723217497714 

0 ( 77 ) = [-59.560673288998(772 - 1) - 9.116379402803(?73 - 1)- 
- 142 . 753455187914(774 - 1) -f 44.843714524611(775 - l)-h 
-f 168.833421156648(77® - 1) - 18.843109698135(771’ - 1)- 
-183.399486372212(77* - 1) -f 2.145923952284(77® - 1)+ 

-f 122.114218185962(7710 - 1) - 36.681425691368(7712 - 1)- 
-0.061343981363(?7i5 - 1) -f 2.491809125293(7714 - 1)] - IO -12 

Case 5.2 For a = ■^, H = 250, the parameters Ci are: 


Cl = -0.011565849071, C 2 = 0.707159448177, C 3 = -290.324617452708, 
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Ci = 435.512207098156, C 5 = -98.780455208880, 

Ce = -0.371954495218, C 7 = -82.314218258861, 

Cg = -37.285501459040- 10■^^ Cg = -15.253404678160 • 10"^^ 
Cio = 12.572809436948- lO'^^ Cn = -22.083271773195 - 10-l^ 
C 12 = 1.090339198400- lO’^^ C 13 = 14.241169061094-10-^2 


and therefore the approximate solutions, (1551) and (H51) may be written as: 


F[ 7 ^) = 1 - 1.638305627622772+ 1.20600966962077^ + 0.043648374556775- 
-1.078984595104776 + 0.15629614192977^ + 0.738925954400778- (52) 

-0.54997260057077^ + 0.122598968736?7io - 0.00021628594677“ 

e{r]) = [-1249.857282929460(772 - 1) - 9.116379418809(773 - 1)- 
-956.378989402938(77^ - 1) + 1351.297888810558(775 - 1)- 
-1018.834483075130(776 - 1) - 646.214253910126(77^ - 1)+ 

+ 1232 . 314945834321(778 - 1) + 129 . 244515099455(779 - 1)- 
-589.041009734213(771° - 1) + 116.176788092689(7712 - l)- 
- 0 . 598803450274(7713 - 1) - 6.703807283538(?7“ - 1)] - IO -12 

Case 5.3 For a = ^, 77 = 500 we obtain: 

+( 77 ) = 1 - 1.1724778890772 + 0.4686048815774 - 0 . 098961423266775 - 
- 0.3550702077776 + 0.78891879354377’’ - 1.8412450049?78+ (54) 

+2.16671419141677° - 1.2119558765771° + 0.25547253506077“ 

%) = [-3025.399926380127(772 - 1) - 9.116379434815(773 - 1)- 
-3566.696711595512(774 - 1) + 4417.102370290691(775 - 1)- 
-2124.393833887592(77° - 1) - 2172.379321553169(77’’ - 1)+ 
+2808.629798261231(778 - 1) + 471.928890384626(77° - 1)- 
-1267.595946519067(771° - 1) + 212.551207395542(7712 - 1)- 
- 1 . 000321320275(7713 - 1) - 11.126037049812(7714 - 1)] - IO -12 
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Case 5.4 For a = ^, H = 1000, it holds that: 


F{r]) = 1 - 0.6223664982r;2 - 0.1660109481??"* - 0.2259232591775+ 
+0.415889430872?76 - 0.57584230277?^ + 0.19874879405077S+ (56) 

+0.0038449003r?9 + 0.0293112856757?5o - 0.05765140267?“ 

%) = [-8164.566371333924(7?2 - 1) - 9.116379466826(7?3 - 1)- 
-12276.720695370957(7?4 - 1) + 13537.207073598438(7?5 - 1)- 
-5217.969886112186(7?6 - 1) - 6768.312414390118(7?^ - 1)+ 

+8081.157802168993(7?® - 1) + 1537.402367705613(7?9 - 1)- 
-3743.086450834657(7?io - 1) + 613.584536626787(7?i2 - l)- 
-1.303403201992(7?i3 - 1) - 31.505440255694(7?“ - 1)] • 10"“ 

175 Case 5.5 For a = ^, = 0, we obtain 

F{r]) = 1 - 1.7695466647r?2 + 1.27541408547?'* + 0.11030235977?5- 
-1.15502567367?® + 0.44340518247?^ + 0.33585643997?®- (58) 

-0.33251921857?® + 0.10455568497?“ - 0.01244219577?“ 

0{r,) = [-111.208165621670(7?2 - 1) - 6.895054111348(7?® - 1)- 
-219.101666428999(7?4 - 1) + 112.000412353392(7?® - 1)+ 
+120.895605108821(7?® - 1) - 47.637714830801(7?^ - 1)- 
-104.355375516421(7?® - 1) + 6.984672282609(7?® - 1)+ 
+75.737792150272(7?“ - 1) - 24.982664253556(7?“ - 1)- 
-0.097447216002(7?*® - 1) + 1.768576986049(7?“ - 1)] • lO-*^ 

Case 5.6 If a = ^, 77 = 250 then 

36 ’ 

F{r]) = 1 - 1.51110142767?® + 0.86128769967?* + 0.01284071907?®- 
-0.56687595097?® + 0.03675984587?*" + 0.31986678997?®- (60) 

-0.16562533477?® + 0.00199 779687?*® + 0.01084986207?** 


13 
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%) = [-5094.867741624686(r;2 - 1) - 6.895054147361(r73 - 1)- 
-4500.773385702808(77^ - 1) + 6249.12145418685(77’5 - 1)- 
-4283.444052831423(77® - 1) - 2834.31757828637(77^ - 1)+ 
+5102.615992457279(77® - 1) + 528.701361603855(77® - 1)- 
-2382.96513866767(7710 - 1) + 460.279590044557(771® - 1)- 
-2.850139405026(7713 - 1) - 26.598784857626(7714 - 1)] • IQ-i® 
Case 5.7 For a = H = 500 the approximate solutions are 


F{ri) = 1 - 1.294221472477® + 0.5334630869774 + 0.001970850477®- 
-0.333060240477® - 0.0197377497771’ + 0.210870962777®- (62) 

-0.115109383677® + O.OI 4765 OO 2 O 77 IO + O.OOIO 58944 O 77 I 1 

0 ( 77 ) = [-10479.441564298402(77® - 1) - 6.895054183374(r7® - 1)- 
-9629.602717192312(774 - 1) + 12532.057010176612(77® - 1)- 
-7708.77471180501(77® - 1) - 5687.792565325205(771" - 1)+ 
+9313.682156269786(77® - 1) + 1063.274966246270(77® - 1)- 
-4206.734533154069(774® - 1) + 761.164805945711(774® - 1)- 
-5.659649234867(774® - 1) - 42.382058179415(7744 - l)] . 10-42 
Case 5.8 li a = ^. H = 1000 then 

36 ’ 


F{r^) = 1 - 0.958190058377® + 0.0913634769774 + 0.003007328577®- 
-0.161106798577® + O.O 66848622677 I’ - 0.069697375977®+ (64) 

+0.045437426277® - 0.0111066789774® - 0.0065559425774® 

0 ( 77 ) = [-22457.262743604428(77® - 1) - 6.895054255400(77® - 1)- 
22247 . 17539598977(774 - 1) + 26605.166448464304(77® - 1)- 
-14220.534699997283(77® - 1) - 12131.493178637833(771’- 1)+ 

(65) 

+ 17921.197074275697(77® - 1) + 2301.819882082563(77® - 1)- 
-7836.564793794578(77®® - 1) + 1301.948445444968(77®® - 1)- 
-11.186604132870(77®® - 1) - 68.576887703688(7744 - 1)] ■ lO"®® 
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From the Tables 1-16 it is obvious that the second-order approximate so¬ 
lutions obtained by OHPM are of a high accuracy in comparison with homo- 
topy perturbation method and with numerical solution obtained by means of 
190 a fourth-order Runge-Kutta method in combination with the shooting method 
using Wolfram Mathematica 6.0 software. 

In Figs 2 and 3 are presented the effect of the Hartmann number on the 
velocity profile for Re = 50 and a = ^ and a = ^ respectively. It is observe 
that velocity increases with increasing of the Hartmann number for any value of 
195 a. The same effect of Hartmann number on the thermal profile are presented in 
Figs 4 and 5 for a = ^ and a = ^ respectively. In this case, the temperature 
decreases with increasing of the Hartmann number in the both cases. The effect 
of the half angle a on the velocity profile is presented in Figs 6-9. With an 
increasing value of a, velocity decreases for H = 0 and H = 250, but increases 
200 for H = 500 and H = 1000. From the Figs 11-13, it is interesting to remark 
that the temperature increases whereas the half angle a increases. In all cases, 
the maximum of temperature occurs near the walls for H = 0 and precisely in 
wall for H ^ 0, while the minimum occurs near the channel axis. 


205 6. Conclusions 

In present paper, the Optimal Homotopy Perturbation Method (OHPM) 
is employed to propose a new analytic approximate solutions for MHD Jeffery- 
Hamel flow with heat transfer problem. Our procedure does not need restrictive 
hypotheses, is very rapid convergent after only two iterations and the conver- 
210 gence of the solutions is ensured in a rigorous way. The cornerstone of the 
validity and flexibility of our procedure is the choice of the linear operator and 
the optimal auxiliary functions which contribute to very accurate solutions. 
The parameters which are involved in the composition of the optimal auxiliary 
functions are optimally identified via various methods in a rigorous way. Our 
215 technique is very effective, explicit, easy to apply which proves that this method 
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is very efficient in practice. 
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Table 1. Comparison between the HPM results [9|, OHPM results (15(1)) and 
numerical results for the velocity F{ri) for a = ^ and H = 0 


r) 

■P’hpm 12] 

F 

^ numeric 

^OHPM" US 

relative error = 

1 ^numeric “ ^OHPmI 

0 

1 

1 

1 

0 

0.1 

0.9770711 

0.9771426047 

0.9771444959 

1.8 '10“® 

0.2 

0.9112020 

0.9114792278 

0.9114802054 

CD 

o 

1 

0.3 

0.8104115 

0.8110052403 

0.8110054657 

2.2 ■10"’’ 

0.4 

0.6859230 

0.6869148220 

0.6869163034 

1.4 -10“® 

0.5 

0.5498427 

0.5512883212 

0.5512895302 

1.2 -10“® 

0.6 

0.4131698 

0.4151088947 

0.4151091740 

2.7 ■10“’' 

0.7 

0.2846024 

0.2870546320 

0.2870554998 

8.6 -lO”’ 

0.8 

0.1702791 

0.1731221390 

0.1731231521 

1.01 -10“® 

0.9 

0.0744232 

0.0768715756 

0.0768721058 

5.3 -10“’ 

1 

0 

0 

0 

0 


Table 2. Comparison between OHPM results (l5T]l and numerical results for 
the temperature 9{ri) for a = ^ and H = 0 


r) 

ft 

'^numeric 

®OHPM mJ 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-9.134405103300 10”’^ 

-9.134559900000 -10“’® 

1.5 -10“’® 

0.1 

-8.553981690585 10“’^ 

-8.561731325450 -10“’® 

7.7 -10“’® 

0.2 

-7.029011445513 -10“’^ 

-7.029011445531 -10“’® 

1.7 10“®® 

0.3 

-4.968059593695 10“’^ 

-4.959904647105 -10“’® 

8.1 -10“’® 

0.4 

-2.830354806908 10“’^ 

-2.830354806921 -10“’® 

1.3 10“®® 

0.5 

-1.003581673783 10“’^ 

-1.015625108634 -10“’® 

1.2 -10“’^ 

0.6 

2.755715234490 10“’® 

2.755715234410 -10“’® 

7.9 10“®^ 

0.7 

9.398183452650 10“’® 

9.682947296376 -10“’® 

2.8 -10“’^ 

0.8 

1.062337244632 10“’® 

1.062337244629 -10“’® 

3.3 10“®^ 

0.9 

7.684282928672 -10“’® 

6.524090716991 -10“’® 

1.1 -10“’® 

1 

0 

0 

0 


19 

































Table 3. Comparison between the HPM results [9|, OHPM results (15^ and 
numerical results for the velocity F{ri) for a = ^ and H = 250 


r) 

■P’hpm 12] 

F 

^ numeric 

^OHPM" HU 

relative error = 

1 ^numeric “ ^OHPmI 

0 

1 

1 

1 

0 

0.1 

0.9837340 

0.9837367791 

0.9837369246 

o 

1 

0.2 

0.9363350 

0.9363459948 

0.9363459260 

6.8 -10“® 

0.3 

0.8616894 

0.8617133581 

0.8617132189 

1.3 •10"’’ 

0.4 

0.7653405 

0.7653814753 

0.7653813980 

O 

1 

00 

0.5 

0.6533961 

0.6534573087 

0.6534570226 

2.8 •10“’' 

0.6 

0.5314621 

0.5315466609 

0.5315462975 

3.6 -10“’ 

0.7 

0.4038130 

0.4039227354 

0.4039224205 

CO 

O 

1 

0.8 

0.2728708 

0.2729980317 

0.2729975130 

5.1 -10“’ 

0.9 

0.1389433 

0.1390416079 

0.1390411070 

5.0 -10“’ 

1 

0 

0 

0 

0 


Table 4. Comparison between OHPM results (l53ll and numerical results for 
the temperature 9{ri) for a = ^ and H = 250 


r) 

ft 

'^numeric 

SqhPM from Eq. 1531 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-8.520036945014 -10“’“ 

-8.520036944914 -10“’° 

9.9 ■10“’'’ 

0.1 

-8.394493067234 10“’° 

-8.395974340239 -10“’° 

1.4 10“’° 

0.2 

-8.032504411694 10“’° 

-8.032504411616 -10“’° 

7.8 ■10“°’ 

0.3 

-7.450795118612 -10“’° 

-7.450297754118 -10“’° 

4.9 10“’^ 

0.4 

-6.677062404515 10“’° 

-6.677062404494 -10“’° 

2.1 10“°’ 

0.5 

-5.746573136109 10“’° 

-5.746735798562 -10“’° 

1.6 ■10“’^ 

0.6 

-4.698328124476 -10“’° 

-4.698328124489 -10“’° 

1.3 10“°’ 

0.7 

-3.570882806332 10“’° 

-3.571240799348 -10“’° 

3.5 10“’^ 

0.8 

-2.398013591946 10“’° 

-2.398013591937 -10“’° 

9.8 10“°° 

0.9 

-1.206024561793 -10“’° 

-1.200670094842 -10“’° 

5.3 10“’° 

1 

0 

0 

0 
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Table 5. Comparison between the HPM results [9|, OHPM results (15^ and 
numerical results for the velocity F{ri) for a = ^ and H = 500 


r) 

■P’hpm 12] 

F 

^ numeric 

■P’OHPM- Eq- mi 

relative error = 

1 ^numeric “ ^OHPmI 

0 

1 

1 

1 

0 

0.1 

0.9883497 

0.9883196668 

0.9883207994 

1.1 -10“® 

0.2 

0.9537955 

0.9537952479 

0.9538026351 

7.3 -10“® 

0.3 

0.8978515 

0.8978510438 

0.8978610430 

9.9 -10"® 

0.4 

0.8224699 

0.8224690439 

0.8224696000 

5.5 •10"’’ 

0.5 

0.7296817 

0.7296803425 

0.7296719365 

8.4 -10“® 

0.6 

0.6209748 

0.6209728597 

0.6209706881 

2.1 -10“® 

0.7 

0.4966644 

0.4966617866 

0.4966700156 

8.2 -10"® 

0.8 

0.3552115 

0.3552079113 

0.3552097261 

1.8 '10“® 

0.9 

0.1923821 

0.1923775215 

0.1923683642 

9.1 -10“® 

1 

0 

0 

0 

0 


Table 6. Comparison between OHPM results (l55]l and numerical results for 
the temperature 9{ri) for a = ^ and H = 500 


r) 

ft 

'^numeric 

SqhPM from Eq. 1551 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-1.783303641361 10“® 

-1.783303641351 -10“® 

9.9 10“®^ 

0.1 

-1.753016328789 -10“® 

-1.753373570315 -10“® 

3.5 10“^® 

0.2 

-1.666810282194 10“® 

-1.666810282186 -10“® 

7.8 -10“®® 

0.3 

-1.531366088677 -10“® 

-1.531258066872 -10“® 

1.08 -10“^® 

0.4 

-1.356325163065 10“® 

-1.356325163062 -10“® 

2.1 10“®® 

0.5 

-1.152391353408 10“® 

-1.152436518900 -10“® 

4.5 -10“®^ 

0.6 

-9.301091200084 ■10“'^° 

-9.301091200095 -10“^® 

1.1 -10“®® 

0.7 

-6.989817680976 -10“^® 

-6.987973346248 -10“^® 

1 

o 

oo 

0.8 

-4.652100850091 10“^® 

-4.652100850080 -10“^® 

1.1 -10“®® 

0.9 

-2.319273360260 ■10“'^® 

-2.321569863571 -10“^® 

2.2 10“®® 

1 

0 

0 

0 
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Table 7. Comparison between the HPM results [9|, OHPM results (155)) and 
numerical results for the velocity F{ri) for a = ^ and H = 1000 


r) 

■P’hpm 12] 

F 

^ numeric 

^OHPM" EH 

relative error = 

1 ^numeric “ ^OHPmI 

0 

1 

1 

1 

0 

0.1 

0.9937607 

0.9937611413 

0.9937578349 

3.3 '10“® 

0.2 

0.9747886 

0.9747905377 

0.9747871858 

3.3 '10“® 

0.3 

0.9422794 

0.9422838643 

0.9422837661 

bo 

o 

1 

00 

0.4 

0.8947431 

0.8947513541 

0.8947499642 

1.3 -lO”® 

0.5 

0.8297471 

0.8297605618 

0.8297564399 

4.1 -10“® 

0.6 

0.7434756 

0.7434959928 

0.7434941281 

1.8 '10“® 

0.7 

0.6299866 

0.6300164190 

0.6300167936 

CO 

o 

1 

0.8 

0.4799338 

0.4799757286 

0.4799724842 

3.2 -10“® 

0.9 

0.2782789 

0.2783297889 

0.2783280591 

1.7 -10“® 

1 

0 

0 

0 

0 


Table 8. Comparison between OHPM results (I57ll and numerical results for 
the temperature 9{ri) for a = ^ and H = 1000 


r) 

ft 

'^numeric 

®OHPM 113 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-3.885903481818 -10“® 

-3.885903481801 -10“® 

o 

1 

to 

o 

0.1 

-3.804390715058 10“® 

-3.805365047292 -10“® 

9.7 10“'® 

0.2 

-3.575104095175 -10“® 

-3.575104095175 -10“® 

5.7 10“°® 

0.3 

-3.222962704714 -10“® 

-3.222630117770 -10“® 

3.3 10“'® 

0.4 

-2.782964569511 -10“® 

-2.782964569511 -10“® 

4.4 10“®® 

0.5 

-2.293213461698 10“® 

-2.293504856351 -10“® 

2.9 10“'® 

0.6 

-1.790055410191 10“® 

-1.790055410191 -10“® 

3.5 10“®® 

0.7 

-1.302936492610 10“® 

-1.301687731203 -10“® 

1.2 -10“'® 

0.8 

-8.442618285524 ■10“'-® 

-8.442618285525 '10“'® 

1.4 10“®® 

0.9 

-4.107432220830 10“'° 

-4.158108458702 '10“'° 

5.06 -10“'® 

1 

0 

0 

0 
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Table 9. Comparison between the OHPM results (1551) and numerical results 
for the velocity F{ri) for a = ^ and H = 0 


r) 

F 

^ numeric 

^OHPM" GU 

relative error = 

(^numeric “ -^OHPmI 

0 

1 

1 

0 

0.1 

0.9824312364 

0.9824320701 

8.3 •10“’' 

0.2 

0.9312259577 

0.9312265466 

bo 

o 

1 

0.3 

0.8506106161 

0.8506107339 

o 

1 

0.4 

0.7467908018 

0.7467915008 

6.9 ■10“’’ 

0.5 

0.6269481682 

0.6269490072 

8.3 -10“’ 

0.6 

0.4982344464 

0.4982347463 

2.9 -10“’ 

0.7 

0.3669663386 

0.3669668037 

4.6 -lO”’ 

0.8 

0.2381237463 

0.2381245686 

8.2 -10“’ 

0.9 

0.1151519312 

0.1151523953 

4.6 -10“’ 

1 

0 

0 

0 


Table 10. Comparison between OHPM results (1591) and numerical results for 
the temperature 9{ri) for a = -^ and H = 0 


r) 

ft 

'^numeric 

SqhPM from Eq. 1591 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-2.552530215568 10““ 

-2.552530214568 -10““ 

9.9 ■10“’’’ 

0.1 

-2.442980325655 -10““ 

-2.444079060328 -10““ 

1.09 -10““ 

0.2 

-2.143998685224 -10““ 

-2.143998685224 -10““ 

1.9 ■10“’’® 

0.3 

-1.714518106214 -10““ 

-1.713385692946 -10““ 

1 

o 

0.4 

-1.227729149275 -10““ 

-1.227729149275 -10““ 

1.3 ■10“’'® 

0.5 

-7.579058130819 -10““ 

-7.588196380463 -10““ 

9.1 10““ 

0.6 

-3.636709774083 10““ 

-3.636709774083 -10““ 

1.6 ■10“’'® 

0.7 

-8.316369531201 10““ 

-8.099883012752 -10““ 

2.1 ■10““ 

0.8 

6.930669117891 10““ 

6.930669117890 -10““ 

7.8 ■10“’'® 

0.9 

9.721365514984 10““ 

8.898023093848 -10““ 

00 

to 

o 

1 

1 

0 

0 

0 
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Table 11. Comparison between the OHPM results (1601) and numerical results 
for the velocity F{r]) for a = ^ and H = 250 


r) 

F 

^ numeric 


relative error = 

(^numeric “ -^OHPmI 

0 

1 

1 

0 

0.1 

0.9849746347 

0.9849746827 

bo 

o 

1 

00 

0.2 

0.9409030596 

0.9409030371 

2.2 -10“® 

0.3 

0.8706210985 

0.8706210491 

4.9 -10“® 

0.4 

0.7783094304 

0.7783094038 

2.6 -lO”* 

0.5 

0.6688194854 

0.6688193877 

o 

1 

00 

0.6 

0.5469607278 

0.5469605971 

1.3 •10“’’ 

0.7 

0.4168757453 

0.4168756336 

1 

O 

0.8 

0.2815737222 

0.2815735380 

bo 

o 

1 

0.9 

0.1426291045 

0.1426289244 

bo 

o 

1 

1 

0 

0 

0 


Table 12. Comparison between OHPM results (IHTl) and numerical results for 
the temperature 9{rf) for a = -^ and H = 250 


r) 

ft 

'^numeric 

®OHPM EJ 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-3.397742006028 10“® 

-3.397742006018 -10“® 

9.9 10“®^ 

0.1 

-3.346638363267 10“® 

-3.347192325339 -10“® 

5.5 -10“^® 

0.2 

-3.199501303768 10“® 

-3.199501303768 -10“® 

1.9 10“®® 

0.3 

-2.964278501171 -10“® 

-2.964072110215 -10“® 

2.06 -10“^® 

0.4 

-2.653181572584 -10“® 

-2.653181572584 -10“® 

6.7 10“®® 

0.5 

-2.281170448213 -10“® 

-2.281224162207 -10“® 

5.3 ■10“®"‘ 

0.6 

-1.864037686164 10“® 

-1.864037686164 -10“® 

1.4 10“®® 

0.7 

-1.416766717489 10“® 

-1.416988166987 -10“® 

2.2 10“®® 

0.8 

-9.521674889800 ■10“'^® 

-9.521674889800 -10“^® 

2.8 10“®^ 

0.9 

-4.798620033623 ■10“'^° 

-4.772452449206 '10“^° 

2.6 10“®® 

1 

0 

0 

0 
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Table 13. Comparison between the OHPM results (1621) and numerical results 
for the velocity F{ri) for a = ^ and H = 500 


r) 

F 

^ numeric 


relative error = 

(^numeric “ -^OHPmI 

0 

1 

1 

0 

0.1 

0.9871408049 

0.9871108182 

1.3 -10“® 

0.2 

0.9490642316 

0.9490642266 

5.0 -lO-® 

0.3 

0.8876104540 

0.8876104486 

5.3 ■10"'* 

0.4 

0.8053144512 

0.8053144616 

1.03 10“® 

0.5 

0.7051032297 

0.7051032241 

5.5 -lO-® 

0.6 

0.5897534597 

0.5897534552 

4.4 -lO-® 

0.7 

0.4613867287 

0.4613867398 

O 

1 

00 

0.8 

0.3210064088 

0.3210063961 

1.2 -10“® 

0.9 

0.1680649652 

0.1680649814 

1.6 -10“® 

1 

0 

0 

0 


Table 14. Comparison between OHPM results (l63l) and numerical results for 
the temperature 9{ri) for a = -^ and H = 500 


r) 

ft 

'^numeric 

SqhPM from Eq. 1631 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-6.861779213872 10“® 

-6.861779213862 -10“® 

9.9 10“®^ 

0.1 

-6.756689795935 10“® 

-6.757837516749 -10“® 

1.1 -10“^® 

0.2 

-6.454596023395 10“® 

-6.454596023395 -10“® 

1.4 10“®® 

0.3 

-5.973076729675 10“® 

-5.972618579646 -10“® 

4.5 10“®® 

0.4 

-5.338639377096 10“® 

-5.338639377096 -10“® 

5.7 -10“®^ 

0.5 

-4.583265280198 10“® 

-4.583356882876 -10“® 

9.1 -10“®^ 

0.6 

-3.739726496008 10“® 

-3.739726496008 -10“® 

6.6 10“®^ 

0.7 

-2.838920429095 10“® 

-2.839180849233 -10“® 

2.6 10“®® 

0.8 

-1.906149341683 10“® 

-1.906149341683 -10“® 

5.7 -10“®^ 

0.9 

-9.597066242185 ■10“'^° 

-9.553857275377 -10“^° 

4.3 10“®® 

1 

0 

0 

0 
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Table 15. Comparison between the OHPM results (1641) and numerical results 
for the velocity F{r]) for a = ^ and H = 1000 


r) 

F 

^ numeric 


relative error = 

(^numeric “ ^OHPmI 

0 

1 

1 

0 

0.1 

0.9904272110 

0.9904271107 

1.003 ■10“'^ 

0.2 

0.9618100677 

0.9618099299 

1.3 ■10“’’ 

0.3 

0.9144036356 

0.9144036639 

2.8 -10"® 

0.4 

0.8484736891 

0.8484737167 

2.7 -lO”* 

0.5 

0.7640642211 

0.7640640849 

1.3 ■10“’’ 

0.6 

0.6606772356 

0.6606771839 

5.1 -10“® 

0.7 

0.5368520578 

0.5368521821 

to 

o 

1 

0.8 

0.3896018441 

0.3896017829 

6.1 -10“® 

0.9 

0.2136111325 

0.2136111294 

3.1 -lO-® 

1 

0 

0 

0 


Table 16. Comparison between OHPM results (1651) and numerical results for 
the temperature 9{rf) for a = -^ and H = 1000 


r) 

a 

'^numeric 

SqhPM fro™ Eq. 1651 

relative error = 

1 ^numeric ~ ^OHPmI 

0 

-1.406496797937 10“® 

-1.406496797936 -10“® 

9.9 -10“®^ 

0.1 

-1.383991714524 10“® 

-1.384237616580 -10“® 

2.4 -10“^® 

0.2 

-1.319483359510 10“® 

-1.319483359510 -10“® 

8.2 -10“®^ 

0.3 

-1.217244871718 -10“® 

-1.217139646331 -10“® 

1.05 -10“^® 

0.4 

-1.083591345359 10“® 

-1.083591345359 -10“® 

1.1 -10“®® 

0.5 

-9.260353523538 10“® 

-9.260364501045 -10“® 

1.09 -10“®^ 

0.6 

-7.519751014022 10“® 

-7.519751014022 -10“® 

4.9 -10“®^ 

0.7 

-5.683228753360 10“® 

-5.683311365697 -10“® 

8.2 -10“®^ 

0.8 

-3.802282135607 -10“® 

-3.802282135607 -10“® 

3.3 -10“®^ 

0.9 

-1.908111552553 10“® 

-1.902733734257 -10“® 

5.3 -10“®® 

1 

0 

0 

0 


26 

































Fig 1. Geometry of the MHD Jeffery-Hamel problem. 



Fig. 2 Effect of the Hartmann Fig. 3 Effect of the Hartmann 

number on the velocity profile for number on the velocity profile for 
a = 7r/24, Re = 50: a = 7r/36, Re = 50: 

-numerical solution, -numerical solution, 



Fig. 4 Effect of the Hartmann Fig. 5 Effect of the Hartmann 

number on the thermal prohle for number on the thermal prohle for 

a = 7r/24, Re = 50: a = 7r/36, Re = 50: 

-numerical solution, -numerical solution, 

. OHPM solution . OHPM solution 
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F(j;) 



Fig. 6 Velocity profile for a = 7r/36, 
and a = 7r/24, H = 0, Re = 50, 
Eqs. (1501) and (1551) : 

-numerical solution, 

. OHPM solution 

Hi) 



Fig. 8 Velocity profile for a = 7r/36, 
and a = 7r/24, H = 500, Re = 50, 
Eqs. (1501) and (1001) : 

-numerical solution, 

. OHPM solution 


Hi) 



Fig. 7 Velocity profile for a = 7r/36, 
and a = 7r/24, H = 250, Re = 50, 
Eqs. (150|) and (1001) : 

-numerical solution, 

. OHPM solution 


Hi) 



Fig. 9 Velocity profile for a = 7r/36, 
and a = 7r/24, H = 1000, Re = 50, 
Eqs. (1501) and (1001) : 

-numerical solution, 

. OHPM solution 
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~eiri) 



Fig. 10 Thermal profile for a = 7r/36, 
and a = 7r/24, H = 0, Re = 50, 
Eqs. (l5T|l and (|5^ : 

-numerical solution, 

. OHPM solution 

e(n) 


-1.x1 

-2.x1 

-3.x1 

-4.x1 

-5.x1 

-6.x1 

-7.x1 






Fig. 11 Thermal profile for a = 7r/36, 
and a = 7r/24, H = 250, Re = 50, 
Eqs. (|53|l and (IHTI) : 

-numerical solution. 


. OHPM solution 

Sin) 



Fig. 12 Thermal profile for a = 7r/36, 
and a = 7r/24, H = 500, Re = 50, 
Eqs. (l55l) and ((63l) : 

-numerical solution, 

. OHPM solution 


Fig. 13 Thermal profile for a = 7r/36, 
and a = 7r/24, H = 1000, Re = 50, 
Eqs. ((57l) and (l65|) : 

-numerical solution, 

. OHPM solution 
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